Wannier states and Bose-Hubbard parameters for 2D optical lattices 
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We consider the physical implementation of a 2D optical lattice with schemes involving 3 and 4 light fields. 
We illustrate the wide range of geometries available to the 3 beam lattice, and compare the general potential 
properties of the two lattice schemes. Numerically calculating the band structure we obtain the Wannier states 
and evaluate the parameters of the Bose-Hubbard models relevant to these lattices. Using these results we 
demonstrate lattices that realize Bose-Hubbard models with 2, 4, or 6 nearest neighbors, and quantify the extent 
that these different lattices effect the superfluid to Mott-insulator transition. 
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I. INTRODUCTION 

There has been considerable recent interest in the study of 
ultra-cold bosonic atoms in optical lattices. The many favor- 
able attributes of optical lattices, such as the absence of de- 
fects, low noise level, and high degree of experimental con- 
trol, are ideal for precise quantum manipulation. As such, 
this system features in theoretical proposals for quantum com- 
puting (e.g. see [1, 2]), simulating many-body systems [3], 
and studying quantum-phase transitions [4]. Beyond bosons, 
work by Hofstetter et ah [5] suggests using an optical lat- 
tice to facilitate the superfluid transition in a degenerate gas 
of fermionic atoms [5]. 

Experimental studies with optical lattices have demon- 
strated many impressive results, including number- squeezing 
an initially coherent Bose-field using a ID lattice [6]; the su- 
perfluid to Mott-insulator transition in a 3D lattice [7]; the 
collapse and revival of coherence in a matter- wave field [8]; 
controlled collisions using state selective lattices [9]. All of 
these experiments are carried out at sufficiently low tempera- 
tures and in deep enough lattices that the constituent bosons 
occupy only the lowest band, and are well described by a tight- 
binding Bose-Hubbard model [10]. In this model the Bose- 
field is decomposed into a set of localized Wannier states, one 
at each site of the lattice, and the Hamiltonian is characterized 
by a few parameters: (i) the number of nearest neighbors sur- 
rounding each lattice site, (ii) the tunneling strength between 
nearest neighbors, and (iii) the strength of the interaction be- 
tween particles at the same site. 

The ground state of the Bose-Hubbard Hamiltonian exhibits 
a quantum phase transition between superfluid and Mott- 
insulating states, that depends on the relative size of the afore- 
mentioned parameters. The transition to an insulating state 
has been suggested as a method for preparing a fiducial state 
of precisely one atom per site, suitable for quantum informa- 
tion processing. Another important application of optical lat- 



*Electronic address: bblakie@physics . otago . ac . nz 



tices is to produce tightly confining potentials suitable for real- 
izing low dimensional gases. For instance, in [1 1] a ID lattice 
was used to investigate Bose-condensation in a quasi-2D gas, 
and experiments with 2D lattices [12, 13] have produced cor- 
related ID gases in the elongated potential tubes at each lattice 
site. For these systems the tunneling parameters and Wannier 
states characterize the time scale over which each potential is 
isolated from neighboring sites, and transverse wavef unction 
for the particles respectively. 

To date most theory, and experiments for quantum- 
degenerate Bose-gases have used cubic lattices, where the op- 
tical potential is produced by orthogonal sets of counter prop- 
agating light fields. For this case the lattice potential is sep- 
arable, and the non-interacting band structure is given by the 
Mathieu spectrum [14]. In this paper we explore the structure 
and properties of a more general class of two-dimensional op- 
tical lattices for regimes relevant to current experiments. We 
compare and contrast the usual square separable lattice (aris- 
ing from two sets of independent laser fields) with a more 
general non- separable lattice produced from three interfering 
light fields. To characterize the appropriate Bose-Hubbard 
models for the lattices considered we numerically evaluate the 
localized Wannier states from band structure calculations. We 
discuss favorable properties of the 3 beam lattice that make 
it more ideal for investigating the strongly correlated regime. 
In particular, we show that certain 3 beam lattice configu- 
rations provide tighter on-site confinement, control over the 
number of nearest neighbors, and significantly reduced tun- 
neling between sites compared to the counter propagating 4 
beam arrangement. These properties of the 3 beam lattice 
lead to the superfluid to Mott-insulator transition occurring 
for much shallower lattices, and would facilitate taking ex- 
periments deeper into the Mott-insulating regime. We also 
show that by using non-counter propagating light fields in the 
4 beam lattice, the depth at which the Mott-transition occurs 
can be significantly reduced, albeit with the requirement of an 
additional light field. 

Conservative two dimensional optical lattices are an im- 
portant tool for the emerging field of quantum atom-optics, 
with numerous practical applications [3, 12, 13, 15]. We see 
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the calculations presented in this paper as serving two impor- 
tant purposes. First, for people studying the idealized Bose- 
Hubbard model, these results will elucidate the typical range 
of Bose-Hubbard parameters accessible, and demonstrate the 
flexibility of the optical lattice system. Second, we believe 
these results will be important in design considerations for fu- 
ture experiments with ultra-cold gases in optical lattices. 



II. FORMALISM 
A. Optical Lattice 



and in the counter propagating case the lattice site separation 
is A/ 2. As illustrated in Fig. 1(a), we have generalized this 
arrangement from counter propagating beams to consider the 
pairs of light fields propagating at symmetrically distributed 
angles with respect to the coordinate axes. As shown below, 
this generalization allows control over the distance between 
lattice sites. 

The lattice potential for this arrangement can be written as 
VWr) = \v cos(bi • r) + ^V cos(b 2 • r), (1) 

where 





bi = ki - k 2 = 2fesin(0/2)x, 
b 2 =k 3 -k 4 = 2fcsin(0/2)y, 



(2) 
(3) 



Figure 1 : Wave vectors of the laser configurations used to generate a 
2D optical lattice, (a) 4 beam lattice formed by two independent sets 
of plane- wave light fields. Different line types are used to differen- 
tiate the two pairs, (b) 3 beam lattice produced by three interfering 
plane-wave fields. The angles are discussed in the text. 

We consider an optical lattice produced from a set 
of monochromatic plane-waves with propagation directions 
characterized by their wave vectors {kj}. To simplify our 
discussion of 2D optical lattices we take these propagation di- 
rections to all lie in a single plane, which we designate as the 
x-y plane. For simplicity we consider all light fields to be 
linearly polarized along z and of equal intensity, with the fo- 
cus of this work on the control of lattice geometry obtained 
by varying the light field propagation directions. The lattice 
potential experienced by the atoms is proportional to the real 
part of the dynamic polarizability a [22], and the light field 
intensity Iq [23]. Below we discuss the two different laser 
configurations we compare in this paper. 

4 Beam Square Lattice In Fig. 1(a) we show the geome- 
try used to realize the 2D optical lattice we refer to as the 4 
beam (4B) square lattice. This type of lattice (and its 3D gen- 
eralization) has been extensively employed in recent experi- 
ments for the case where each pair of lasers is counter prop- 
agating (i.e. with <j> = 180°) [7, 15]. The two pairs of light 
fields are made independent from each other by detuning the 
common frequency in one pair of fields from that of the other 
pair. Typically, a negligible frequency difference compared 
to the optical frequency is required to achieve independence 
(e.g. see [16]), so to a good approximation all wave- vectors 
are the same magnitude |kj| w k = 2tt/X j where A is the 
optical wavelength. The resulting potential is a superposition 
of two perpendicular ID lattices (one arising from each pair), 



are the reciprocal lattice vectors and Vb = — Re{a}/ / 2e c is 
four times the dipole shift associated with a single laser of in- 
tensity io- We will refer to \ Vo \ as the light-shift strength, and 
note that \ Vo\is the potential depth for the ID lattice produced 
by two intersecting light fields of equal intensity Io . In more 
than one spatial dimension the potential saddle point between 
sites is the most essential characterization of the lattice depth, 
however this depends strongly on the lattice geometry and is 
difficult to characterize. To be of most relevance to experi- 
ments, in this paper we will compare different lattices made 
with the same light- shift strength Vb, i.e. produced by sets 
of laser fields with the same intensity. We note that Vq can 
be a positive or negative quantity depending on whether the 
light fields are blue or red detuned from atomic resonance re- 
spectively. The 4 beam lattice is symmetric with respect to Vb 
changing sign (to within an overall translation of the lattice), 
however the next lattice we consider is not. 

The direct lattice vectors which specify the translation vec- 
tors between lattices sites are 



^ csc(0/2)x, 



a 2 = -csc(0/2)y, 



(4) 
(5) 



demonstrating that the lattice site spacing can be varied from 
A/2 — > oo as changes from 180° — > 0. 

3 Beam Lattice The second 2D optical lattice we consider 
consists of 3 laser fields — the minimum number needed to 
generate a 2D lattice. We will refer to this configuration as 
the 3 beam (3B) lattice. The geometric arrangement of the 
light fields is shown in Fig. 1(b), with all three light fields 
taken to have the same frequency. The lattice potential, for 
this arrangement, given by 

VLatt(r) = ^Vocos(bi-r) + ivocos(b 2 -r) (6) 
+ iy cos([bi + b 2 ]-r), 

is a sum of three ID lattices arising from the interference be- 
tween each distinct pair of fields. As for the 4 beam case, 
we have introduced the light- shift strength parameter Vq = 
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— Re{a}2o/2eoc, and we define the reciprocal lattice vectors 
bj below. 

The incident wave- vectors are of the form 



(7) 



and generate a potential that exhibits a broad range of 
crystallographic configurations, dependent on the choice of 
{#1, #2, #3}. The crystallography of this type of lattice has 
been considered by Petsas et al. [17] in a context relevant to 
laser cooling atoms, however we find it useful to cast their ob- 
servations in an analytic framework relevant to conservative 
lattices. To do this we choose to label our wave vectors such 
that 0i < 02 < #3, and place the following restrictions on the 
angles: 



i=0, 
2 < 180°, 



< 

< 3 - 02 < 180°, 



(8) 
(9) 
(10) 



For any set of three coplanar light fields (excluding the case 
where two beams are co-propagating which gives rise to a ID 
lattice) the angles can always be transformed to satisfy re- 
quirements (8)-(10) with an overall rotation the system in the 
x-y plane. 

Analogously to what was done in the four beam lattice, we 
define the reciprocal lattice vectors as 



(11) 



We note that other choices of reciprocal lattice vectors are 
possible. The direct lattice vectors {ai, ^2} can then be de- 
termined from the orthogonality relationship a^ • bj = 2ir5ij 
(e.g. see [18]). For the purposes of discussing direct lattice 
geometry, it is useful to specify the angle </>£> between the di- 
rect lattice vectors and the ratio of their lengths vr>. It can be 
shown that these relate to the incident laser angles as 



td = 



PD 
1*2 I 



180°-^, 
sin(<9 2 /2) 



|ai| sin([0 3 - 2 ]/2)' 



(12) 
(13) 



This expression can be inverted to uniquely determine the 
laser angles under restrictions (8)-(10) for m £ (0, 00), and 
<t>D e (cos-^l/r^lSO ). 

Finally, we note that unlike the 4 beam potential, the 3 beam 
potential (6) is not symmetric with respect to sign change of 
Vo (i.e. change in the sign of the laser detuning from atomic 
resonance), a point we elaborate on in Sec. III. The 3 beam 
lattice only realizes a periodic set of confining wells for the 
case of red detuning, i.e. Vo < 0, and as such we restrict our 
attention to red-detuned lattices in this paper. 



B. Wannier states 

The eigenstates of the periodic single particle Hamiltonian, 

^ = p 2 /2m + F Latt (r), (14) 



are known as Bloch states, which we write as ^ nq (r), with 
respective energy eigenvalue ?kj nq . Here we take the lattice 
potential to be a general periodic function on the Bravais lat- 
tice of points Rj = nij3.i + U2j8L2 {^ij £ such that 
VLatt(r) = V La tt(r + Rj). According to Bloch's theorem, 
these eigenstates can be written as 



VW r ) = ^nq(r)e* q,r , 



(15) 



where u ncL (r) 



(r + a$) is a periodic function, q is 
the quasi-momentum, and n is the band index. Formally, a 
Bravais lattice is of infinite extent, however it is numerically 
convenient for us to consider a finite lattice of N sites with 
periodic boundary conditions. For this situation each energy 
band contains TV eigenstates, and the quasi-momenta can be 
restricted to the first Brillouin zone. 

Bloch states are extended states that span the entire lattice. 
To consider many-body effects arising from the interactions 
between particles it is more convenient to work with a set of 
localized states. Here we consider the Wannier basis which is 
a unitary transformation of the Bloch basis, where the basis 
states are labeled by lattice site position. The Wannier state 
centered at site R^ is defined as 

w n (r-K t ) = ^^ e -^ R ^nq(r), (16) 
* q 

where the summation is taken over all Bloch states of band- 
n. For describing quantum degenerate Bose-gases in optical 
lattices we will only consider Wannier states of the ground 
band. 

Wannier states are not eigenstates of the Hamiltonian and so 
an atom prepared into a Wannier state at a given site will tun- 
nel to other Wannier states over time. This tunneling between 
sites is characterized by the matrix elements of the Hamilto- 
nian between the respective Wannier states at those sites, i.e. 

7Ri-Rj = J drw^r - K^HqWq^ - TLj), (17) 



which reduces to 



-iq-(Rj— Ri) 



(18) 



i.e. the Fourier transform of the Bloch dispersion relation. 
We note that the overlap integral between states at different 
sites (18) only depends on the relative separation between the 
sites, so that characterizing the tunneling properties of a single 
Wannier state completely determines the matrix elements of 
the whole lattice. 

For deep lattices the Wannier state is well localized at its 
central site and 7r rapidly decreases with increasing relative 
separation |R|. This is the well-known tight binding limit, 
where a state at any given site only couples (to any degree 
of significance) to a few neighboring sites. In this limit only 
the most prominent tunneling matrix elements to other sites 
are retained, where these sites are referred to as the nearest 
neighbors. For cases where the nearest neighbor tunneling 
terms are all equal it is conventional to define the tunneling 



4 



parameter or hopping rate as J = — 7rj (e.g. see [4]). We 
will present quantitative results for these tunneling rates and 
demonstrate systems which exhibit different numbers of near- 
est neighbors in Sec. IV. 



C. Bose-Hubbard model 

Here we review the Bose-Hubbard model for interacting 
Bosonic atoms in an optical lattice, closely following the orig- 
inal derivation by Jaksch et al [4]. The usual starting point 
for discussing many-body aspects of ultra-cold atoms is the 
second quantized Hamiltonian 

H = ^r^t( r )|^ + y ext ( r )}^( r ) 

+ ±g j dv^{v)^{v)^{v)^{v), (19) 

where -0(r) is the boson field operator, ^4 x t( r ) describes a 
slowly varying external potential (in addition to the lattice po- 
tential in Hq), g = Anh 2 a s /m is the interaction strength, with 
a s the scattering length and m the atomic mass. 

If a sample of sufficiently cold Bosonic atoms is loaded into 
a deep enough lattice (e.g. from a Bose-Einstein condensate), 
then only states of the lowest vibrational band will be popu- 
lated. This system is conveniently described by expanding the 
field operators in the Wannier basis of the ground band 



^(r) 



w (r - TLi)a,i 



(20) 



where the operator di destroys a particle in the -Wannier 
orbital, and satisfies the usual Bose commutation relations 
[ai^aj} = and [ai,aj] = Sij. In the Wannier representa- 
tion the Hamiltonian (19) reduces to the Bose-Hubbard form 



(21) 



where the first summation in Eq. (19) has been restricted to 
the nearest neighbor sites (i.e. tight binding approximation 
discussed above). We have also introduced the site dependent 
local energy 

Ci = JdTw^(r-K i )V ext (T)w (r--R i ), (22) 
- Kxt(Ri) (23) 
and the on-site interaction term 

U B n=gj dr\w (r)\\ (24) 

The Bose-Hubbard model will be a good description of 
the system if degrees of freedom in the higher bands can 



be neglected. For an energy gap to the first excited band of 
^gap, we find the validity conditions to be /c#T <C E gap and 
uUbh <C i^gap* where n is the average site occupation. Typi- 
cally these requirements are well satisfied in experiments. 

In this paper we calculate 2D Bose-Hubbard model parame- 
ters that can be achieved using a 2D optical lattice. The nature 
of the confinement transverse to the lattice will have consider- 
able influence on the physical properties of the system, how- 
ever this is beyond the scope of this paper. For simplicity we 
assume that the motion transverse to the lattice is frozen out, 
e.g. as could be arranged experimentally by superimposing a 
deep ID lattice to slice the potential tubes of the 2D lattice 
into a set of independent planes. 



D. Superfluid to Mott-insulator transition 

One of the most fascinating properties of the Bose-Hubbard 
Hamiltonian is that it predicts a quantum-phase transition of 
the ground state of the system between superfluid and Mott- 
insulating states [10]. The phase diagram for the system de- 
pends precisely on the Bose-Hubbard parameters and the av- 
erage filling factor of the system. It is convenient to define the 
dimensionless quantity 



zJ ' 



(25) 



embodying the Bose-Hubbard model parameters, where z is 
the number of nearest neighbors. For a homogeneous system 
with unit filling (i.e. an average density of one boson per site), 
a critical value of r] c w 5.8 is predicted by mean-field the- 
ories for the transition [24]. For rj < rj c the system is in a 
superfluid state where the hopping between sites dominates 
over on-site repulsion. In this state, the system exhibits phase 
coherence across the lattice and significant fluctuations in the 
number of bosons per site. For r] > f] c interaction effects dom- 
inate tunneling and the particles become localized — this is 
the Mott-insulator state, characterized by a near definite num- 
ber of particles per site and an absence of phase coherence 
between sites. 

For bosonic atoms in an optical lattice the values of Ubh 
and J can be dynamically modified by changing the lattice 
depth. E.g. the effect of increasing the lattice depth is to con- 
fine the bosons more tightly at each site, increasing Z7bh> and 
reducing the tunneling between sites. Greiner et al [7] have 
experimentally demonstrated a loading procedure that is suffi- 
ciently adiabatic for the system to be reversibly taken through 
the phase transition. 



III. LATTICE GEOMETRY 
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Figure 2: Optical lattice potentials, (a)-(e) Potentials for the 3 beam laser configuration, (f) Potential for the 4 beam counter propagating 
configuration. For these plots Vb = —9Er, other parameters discussed in the text. 



In Fig. 2 we demonstrate a variety of the possible 2D opti- 
cal lattice geometries that can be made with the lattice types 
discussed in Sec. II A. In particular, the 3 beam configuration 
allows the five distinct types of two-dimensional Bravais lat- 
tices to be realized by suitably choosing the incident light field 
angles (Figs. 2(a)-(e)). In Fig. 2(f) we show the counter prop- 
agating 4 beam lattice potential. Note that we have chosen 
to specify light- shift strength in units of photon recoil energy, 
defined as E R = h 2 k 2 /2m. 

In more detail the cases considered in Fig. 2 are (named 
according to the underlying Bravais lattice geometry): 

• (a) [3 Beam] Hexagonal Lattice: This highly symmet- 
ric triangular lattice displays hexad symmetry (symme- 
try under 60° rotations). This lattice naturally arises in 
close packing of spheres - and corresponds to the planes 
of three dimensional hexagonal close-packed and face- 
centered cubic lattices [18]. The direct lattice vectors 
have length ratio ro = 1 and intersect at an angle of 
<pD = 60°. The light fields needed to make this lattice 
have planar angles of 2 = 120.0° and <9 3 = 240.0° 
(seeEqs. (12) and (13)). 

• (b) [3 Beam] Square Lattice: This square lattice has an 
underlying Bravais lattice with tetrad symmetry (sym- 
metric under 90° rotations), though the potential does 



not exhibit this symmetry (it is only symmetric under 
180° rotations). The direct lattice vectors have length 
ratio td = 1 and intersect at an angle of 4>d = 90°. 
The light fields needed to make this lattice have planar 
angles of <9 2 = 90.0° and <9 3 = 180.0°. 

• (c) [3 Beam] Rectangular Lattice: The rectangular 
lattice emphasizes the considerable asymmetry between 
directions that can be engineered. For the specific case 
we have chosen to consider the direct lattice vectors 
have length ratio td = 3.5 and intersect at an angle 
of (j)D = 90°. The large lattice vector a 2 arises because 
the light fields with wavevector k 2 and k 3 differ by a 
relatively small angle of 31.9°, i.e. 2 = 148.1° and 
0s = 180.0°. 

• (d) [3 Beam] Oblique Lattice: This is the most general 
planar lattice type of which all other 2D Bravais lattices 
are special cases. The only symmetry that this lattice 
exhibits in general is diad symmetry (symmetric under 
180° rotations). For the specific case considered, the 
direct lattice vectors have length ratio ro = 0.6 and in- 
tersect at an angle of = 85°. The light fields needed 
to make this lattice have planar angles of # 2 = 64.5° 
and 3 = 190.0°. 
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• (e) [3 Beam] Centered Rectangular Lattice: The cen- 
tered rectangular lattice is usually specified using the 
larger (non-primitive) cell indicated in Fig. 2(e), that 
emphasizes the presence of a face centered point. For 
the specific case considered, the direct lattice vectors 
have length ratio m = 1 and intersect at an angle of 
<j)D = 35°. The light fields needed to make this lattice 
have planar angles of 2 = 145.0° and 3 = 290.0°. 

• (f) [4 Beam] Square Lattice: This lattice is formed 
from two sets of counter propagating light fields with a 
A/2 lattice site spacing, and exhibits the highest lattice 
site density (per unit area) of the lattices considered in 
Fig. 2. 

Before considering the Wannier states of these lattices it is 
useful to compare the most significant differences between the 
3 beam and 4 beam lattices. 

Lattice geometry: As demonstrated in Figs. 2(a)-(e), the 
3 beam lattice exhibits a wide range of lattice types, with the 
laser angles calculated using expressions (12) and (13). The 
rigid constraints we placed upon the light field arrangement of 
the 4 beam lattice ensures the resulting lattice is square. Re- 
laxing these constraints by allowing the relative intersection 
angle of each pair to be different, or apply an overall rotation 
to one pair, it is possible to generate rectangular or oblique 
lattices respectively. These generalizations are not considered 
in this work. 

The lattice geometry plays a central role in determining the 
number of nearest neighbors, z, in the derived Bose-Hubbard 
model. Inspection of Fig. 2 suggests that the 3 beam lattice of- 
fers cases that realize 2 (see Fig. 2(c)), 4 (see Fig. 2(b)), and 6 
(see Fig. 2(a)) nearest neighbors. With the aide of band struc- 
ture calculations for the tunneling matrix elements we return 
to this issue in the next section, and verify the actual numbers 
of nearest neighbors. Another important factor relating to the 
geometry is the distance between lattice sites. For the 4 beam 
lattice (Fig. 2(f)) the spacing between sites is the minimum 
value of |aj| = A/ 2. For the 3 beam lattice spacing is gen- 
erally larger, e.g. \a.j\ = 2 A/ 3 for the hexagonal lattice (Fig. 
2(a)), \a.j | = A/ y/2 in the square lattice (Fig. 2(b)). The site 
spacing will influence the tunneling between sites, as this is 
strongly dependent on the barrier width between sites. 

Unit cell potential shape: The single well potential for 
the 3 beam and 4 beam lattices are quite distinctive, as is illus- 
trated for two particular cases in Fig. 3. The most noticeable 
feature distinguishing these cases (and is a general observa- 
tion for the different configurations of 3 beam lattice) is that 



the saddle point of the potential separating a given lattice site 
from the neighboring site is much higher in the 3 beam lattice 
(Fig. 3(a)) than in the 4 beam lattice (Fig. 3(b)). The superior 
confinement of the 3 beam lattice produces much better local- 
ized Wannier states than for a 4 beam lattice with the same lat- 
tice site spacing and field intensity (per incident beam), as we 
demonstrate in the next section. From Fig. 3(a) it is clear why 
the 3 beam lattice is not symmetric with respect to Vb chang- 
ing sign: The subtle crown-like maxima of the 3 beam lattice 
(which become the minima when the light fields are changed 
to blue detuning) are shallow and provide little confinement. 
The dominant feature for the blue detuned 3 beam lattice is 
the periodic array of maxima, referred to as an anti-dot lattice 
(e.g. see [19]). 




Figure 3: Comparison of the individual potential wells of the 3 beam 
and 4 beam lattices. Potential energy contour plot of (a) 3 beam 
hexagonal lattice, and (b) 4 beam square lattice. The potentials are 
displayed on a region of space slightly larger than a primitive unit 
cell. The lattice potentials displayed here correspond to the same 
configurations and parameters used in Figs. 2(a) and (f) respectively. 
An energy offset has been added so that the potential minima are at 
0. 



IV. WANNIER STATES AND BOSE-HUBBARD 
PARAMETERS 



In this section we present numerical calculations of the 
Wannier states and Bose-Hubbard parameters for the differ- 
ent 2D optical lattices considered in this paper. 

To begin, we briefly summarize our numerical procedures. 
We diagonalize the non-interacting lattice Hamiltonian (14) 
represented in a plane-wave expansion to obtain the ground 



band Bloch states and energy eigenvalues for each quasi- 
momentum value within the first Brillouin zone. Ensuring 
the Bloch states are continuous in quasi-momentum, we con- 
struct the Wannier states by numerically summing the Bloch 
states with the appropriate phase factors (see Eq. (16)). Hav- 
ing obtained the Wannier states, the on-site interaction can be 
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Figure 4: Wannier state density distributions and tunneling neighborhoods. The log-densities of the Wannier states in (a)-(f) correspond to the 
similarly labeled potentials displayed in Fig. 2. Potential energy contours are indicated by dotted white lines mark the lattice site locations. 
The white numbers centered upon neighboring sites denote the 'nearness' of that site according to a ranking of the tunneling matrix elements 
(see text). For these plots Vb = —9Er. 



evaluated according to Eq. (24) using numerical integration. 
Finally, the tunneling matrix elements are calculated from the 
ground band energy spectrum using Eq. (18). We note that 
to numerically calculate the tunneling matrix elements cou- 
plings, 7r,_r j , we need to take the lattice size TV to be suffi- 
ciently large to contain the sites and Hj . 

The Wannier state localized at the origin site for each of 
the six lattice potentials of Fig. 2 are shown in Fig. 4. At 
the light- shift strength considered (Vb = —9Er), the Wannier 
states are all well localized upon the central site, and to facil- 
itate observing the small density tails at the neighboring sites 
we have employed a log- scale for the density color-map. The 
Wannier state of the 4 beam square lattice is the least localized 
of the all the states considered, mainly due to poor degree of 
lattice site confinement discussed earlier (see Fig. 3(b)). It 
is interesting to note that the first node in the Wannier states 
occurs near the middle of the neighboring site, as is required 
for states at different sites to be orthogonal. 

The tunneling matrix elements provide a quantitative mea- 
sure of the coupling strength between different Wannier states. 
In order to demonstrate the variable number of nearest neigh- 
bors, and the extent to which further than nearest neighbor 
couplings are negligible, we have listed the tunneling matrix 
elements in table I and have indicated the groups of sites these 



matrix elements correspond to in Fig. 4. To be clear, in Fig. 
4 we have enumerated all sites that have the same tunneling 
matrix element to the central site with the same number, e.g. 
1 indicates 'first nearest neighbor', 2 indicates '2nd nearest 
neighbor', etc. (i.e. larger integers indicate decreasing magni- 
tude of tunneling matrix element). These results indicate that 
at Vb = —9Er the second nearest neighbors have tunneling 
rates that are 2-4 orders of magnitude smaller than those of 
the nearest neighbors. Tight-binding can be expected to pro- 
vide an excellent description of the system in this regime, i.e. 
we are justified in identifying the first column of Table I with 
— J (see Sec. II B), and ignoring all other couplings. Another 
important observation is that the number of nearest neighbors 
for the Wannier states in Figs. 4(a)-(f) (i.e. the number of l's 
in each figure) takes the values 2, 4, or 6. 

We note that the tunneling to the rows of sites vertically sep- 
arated from the site at the origin in Fig. 4(c) is so small that it 
is completely negligible compared to the first three neighbors 
considered. Also, the separable nature of the potential for the 
square 4 beam lattice in Fig. (f) means that tunneling to all 
sites not lying on axis are identically zero. 

As a final comparison between the lattices we consider the 
global variation of the nearest neighbor tunneling and the on- 
site interaction for a broad range of light- shift strengths. We 
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site label: 


1 -sites 


2- sites 


3 -sites 


(a) (3B) Hexagonal 


-2.79 x 1(T 


4 


8.75 x 10- 8 


5.64 x 10~ 8 


(b) (3B) Square 


-1.862 x 10 


-4 


-1.12 x 10" 6 


-8.88 x 10" 8 


(c) (3B) Rectangular 


-2.90 x 10" 


3 


4.54 x 10" 6 


-1.12 x 10" 8 


(d) (3B) Oblique 


-1.83 x 10" 


3 


-6.48 x 10" 6 


-4.28 x 10" 6 


(e) (3B) Centered 


-1.80 x 10" 


3 


-6.33 x 10" 6 


1.84 x 10" 6 


(f) (4B) Square 


-2.42 x 10" 


2 


3.77 x 10~ 4 


-9.28 x 10" 6 



Table I: Values of the first three nearest neighbor couplings for the 
Wannier states considered in Fig. 4 in units of Er. The values in the 
columns labeled "j-sites" correspond to the tunneling matrix element 
between the Wannier state at the origin and the sites enumerated by 
j in Fig. 4. 



have chosen five lattices to make this comparison. The first 
four lattices are the most unique of those considered in Figs. 
2 and 4: (a) the 3 beam hexagonal lattice (6 nearest neigh- 
bors); (b) the 3 beam square lattice (4 nearest neighbors); 
(c) the 3 beam rectangular lattice (2 nearest neighbors); and 
(f) the 4 beam square lattice (4 nearest neighbors). For the 
fifth lattice we introduce a new 4 beam configuration with a 
larger lattice constant of \a. 3 -,\ = 2 A/ 3 (for which the pairs of 
light fields intersect at (/> w 97.2°), chosen to have the same 
distance between sites as the hexagonal lattice. Our motiva- 
tion for considering this lattice is to better assess the extent 
to which wider spacing between sites in the 4 beam lattice 
reduces nearest neighbor tunneling. 





5 10 15 20 25 

l V ol [ E r] 

Figure 5: Bose-Hubbard parameters for a variety of different 2D lattices at various light- shift strengths (a) Nearest neighbor tunneling strength, 
(b) On-site interaction strength, (c) Dimensionless scaling parameter 77 (see Sec. II D). Dashed horizontal line indicates the critical value for 
superfluid to Mott-insulator transition for a homogeneous system with filling factor n = 1. We have taken an effective 2D coupling constant 
of g2D = 0.98 ER/k 2 (k — 2tt/X). Note that the coupling strength g, discussed in Sec. II C, involves the 3D scattering length, so that #2D, 
is obtained by integrating over the degree of freedom transverse to the lattice. The value of g^D we use here assumes the potential transverse 
to the lattice is a tightly confining harmonic trap of frequency cu z = hk 2 /m. We choose physical parameters corresponding to the scattering 
length of 87 Rb and A = 852nm (similar to those of Ref. [7]). 

I 



The results for the five chosen lattice types are shown in ments are displayed. Generally the 3 beam lattices have corn- 
Fig. 5. In Fig. 5(a) the nearest neighbor tunneling matrix ele- 
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parable or considerably smaller tunneling matrix elements 
than the 4 beam lattices at any given light- shift strength. In 
particular, tunneling in the 3 beam square and hexagonal lat- 
tices is more that 3 orders of magnitude smaller than the 
counter propagating 4 beam lattice at the maximum light- 
shift. Of the 4 beam lattices, the non-counter propagating con- 
figuration has significantly smaller tunneling matrix element 
than the counter propagating configuration at each light- shift 
strength. Nevertheless, the non-counter propagating 4 beam 
lattice (with lattice site spacing matching that of the hexagon 
lattice) has much larger tunneling matrix elements than the 
hexagonal lattice due to the residual differences in the poten- 
tials of the individual wells (see Sec. III). 

In Fig. 5(b) we show the on-site interaction strength, Ubyl- 
It is surprising to note that for the hexagonal lattice Ubyl is 
higher than that of the counter propagating 4 beam lattice at 
every value of Vb, indicating that the hexagonal lattice has 
tighter lattice site confinement. A useful measure of the con- 
finement is furnished by the oscillation frequency of the har- 
monic approximation to the lattice site minima. For fixed Vo, 
we find that the oscillation frequency of the hexagonal lattice 
is ~ 6% higher than for the counter-propagating 4 beam lat- 
tice. 

In Fig. 5(c) we plot r\ (see Eq. (25)), which embodies the 
parameters of the Bose-Hubbard Hamiltonian. As discussed 
in Sec. II D, the value of r\ determines whether the equilib- 
rium ground state is a superfluid or a Mott-insulator. For ref- 
erence we have indicated the transition point for filling factor 
n = 1 as a horizontal dashed line in Fig. 5(c). For the lattices 
considered, the Mott-transition (where the curves cross the 
dashed line) varies over a wide range of light- shift strengths. 
In the hexagonal lattice the Mott-transition is at the lowest 
light- shift of | Vb| ~ 3.7 Er, whereas the transition for the 4 
beam counter propagating lattice occurs at the largest value of 
| Vb| ~ 10.2Er. For reference, this 4 beam lattice is the two 
dimensional version of the lattice recently used by Greiner et 
al [7] to experimentally observe the Mott-transition. 

Comparing two 4 beam lattices in Figs. 5(a)-(c) yields the 
following observations. Relative to the |aj| = A/2 lattice, 



the |aj| = 2 A/3 lattice has significantly smaller J, smaller 
/7bh, which combined result in a substantially larger rj value 
at each light-shift strength. In general choosing smaller angles 
between the pairs of light fields <p in the 4 beam lattice (see 
Fig. 1) could be used to maximize rj, and reduce the light- shift 
strength at which the Mott-transition occurs. E.g., for cj> = 
60° the lattice site spacing is A, and Mott-transition occurs 
at | Vo| ~ 2.5Er — lower than any of the 3 beam curves. 
However, reducing the transition light- shift strength is made 
at the expense of Ubji, the tightness of well site confinement, 
and the energy gap to higher bands (hence the validity of the 
Bose-Hubbard description). Balancing these trade offs will be 
an important consideration for many practical applications. 



V. CONCLUSIONS 

In this paper we have investigated different types of 2D op- 
tical lattices that might be used in experiments to produce 
highly correlated states, such as the Mott-insulator state, or 
to engineer ID potentials. Through numerical investigation 
of the band structure and explicit calculation of the Wannier 
states, we have evaluated the Bose-Hubbard parameters over 
a wide regime. We have demonstrated lattices that realize 
Bose-Hubbard models with 2, 4, or 6 nearest neighbors and 
significantly effect the light- shift strength at which superfluid 
to Mott-insulator transition occurs. 
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